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Abstract We theoretically study the nonlinear dynamics of the instability of counter- 
superflow in two miscible Bose-Einstein condensates. The condensates become 
unstable when the relative velocity exceeds a critical value, which is called counter- 
superflow instability. We reveal that the counter- superflow instability can lead to 
quantum turbulence by numerically solving the coupled Gross-Pitaevskii equa- 
tions. The modes amplified by the instability grow into solitons and decay into 
quantized vortices. Eventually, the vortices become tangled and quantum turbu- 
lence of two superfluids. We show that this process may occur in experiments by 
investigating the dynamics in a 2D trapped system. 
PACS numbers: 03.75.Kk,03.75.Mn,05.30.Jp 



1 Introduction 

Quantum turbulence (QT) is currently one of the most important topics in low 
temperature physics ^ QT has been historically studied in superfluid "^He and 
■^He, while QT is now numerically studied in atomic Bose-Einstein condensates 
(BECs)^ and experimentally realized recently-^. Atomic BECs have several ad- 
vantages over superfluid helium. The most important point is that modern optical 
techniques enable us to easily control the condensate and directly visualize quan- 
tized vortices. 

Old studies on QT in superfluid "^He focused primarily on thermal counterflow, 
in which the normal fluid and superfluid flow in opposite directions driven by an 
injected heat current. In this work we consider a similar case in two-component 
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Fig. 1 QT in two-component BECs. This shows the surface of the equal density, n = 0.1. The 
black one and the gray one are the vortices of component 1 and component 2 . 



BECs. Two-component BECs are known to create various exotic structure of 
quantized vortices'^ and cause some characteristic hydrodynamic instability such 
as Kelvin-Helmholtz instability^ and Rayleigh-Taylor instability^. Hence QT and 
the instability towards QT in two-component BECs should propose novel dynami- 
cal phenomena and open the new research field in cold atoms. For example, the en- 
ergy spectra of QT in a single-component BEC is known to obey the Kolmogorov 
law that is the most important statistical law in turbulence but it is not so trivial 
what happens to QT of two-component BECs. 

In this work we study numerically QT in counterflow of two-component mis- 
cible BECs for the first time (Fig. 1). Law et al. showed theoretically that counter- 
superflow is unstable when the relative velocity exceeds some critical value by 
linear analysis^. Very recently Hamner et al. study experimentally and numeri- 
cally shocks and dark-bright solitons in counterflowing BECs^. Our work focuses 
on the drastic state of QT and the scenario of the instability towards QT. 



2 Instability of counter-superflow 

We consider binary mixture of BECs described by the condensate wave functions 
Wj — ^ynje^^j in the mean-field approximation at T = K, where the index j refers 
to each component j (7=1,2). The wave functions are governed by the coupled 
Gross-Pitaevskii (GP) equations 

//^|-^i = -;^v2^i + y(r)^i+gii|^ip^i+gi2|^2|'^i, (1) 

= - V V2 + y (r) 'F2 + g22 1 'ftl V2 + gu I fl I V2, (2) 
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where mj is particle mass associated with the species, gjj is intracomponent inter- 
action and gi2 is intercomponent interaction. Our analysis satisfies the condition 
^11^22 > gu ^^^^ condensates are miscible. Following the parameters in the 
experiments^, we consider the case that particle mass mj and interaction coeffi- 
cient gjj of two condensates are equal, namely mu = m22 = ^ and gu = g22 = ^• 
We can obtain the critical relative velocity in counter- superflow from dynami- 
cal instability. In a stationary state, the wave functions are written as Wj^ = ^^Jr^-e^^y^ , 

where superfluid velocity is Vj — ^ V0j — ^kj. We assume n\ — n2 — n^. The dis- 
turbance 5^- around is represented by 

8Wj = e^^^J-r-^jt) {^^.^^(^-^-^O _ v*^-^(^-^-^^) }. (3) 

Inserting these formulations into the GP-equation yields the Bogoliubov-de Gennes 
equation. Solving the eigenvalue problem, the frequencies CO are found to obey 

{hw-Vc'Pf = £^^2£grP^^\VR'P\^±^{£^^2£gn^)\VR-P\^^4£^g^n^^y^. 

(4) 

Here Vc = ( Vi + V2) /2 is the velocity of the center of mass, Vr = V1 — V2 

is the relative velocity, P — hq is the momentum of excitation, £ = ^ is the 
energy of excitation, and 7 = gu/g- When an imaginary part of a frequency is 
nonzero,/.^, {hco — Vc ■ P)^ < 0, the system is dynamically unstable. Then we 
obtain the condition of instability of counter- superflow is 

4(ii+£r!!_£i!,)<,v4.p)^<4(^i+si! + 5ri!r). (« 

\4m^ m m / \4m^ m m J 

Figure 2 shows this condition with 7 = 0.5. The horizontal axis is the excited 
mode P and the vertical axis is the relative velocity Vr. The gray zone shows that 
counter-superflow is unstable and the P = mode becomes unstable first with 
increasing the relative velocity Vr. Hence the critical relative velocity V^riticai is 

\j m m '' 

The dynamical instability should induce the exchange of momentum between 
the two components to reduce the relative motion of the condensates. Here we may 
define the momentum of the component j as Jj = / drYii^f^^j — ^V^*). In a 

stationary counter-superflow state with = y/^exp(/ ^^^ ^ ), we have /V — 
mrijVj, where we use the volume V of the system. Let us suppose that a single 
unstable mode appears in the counter-superflow state. Then the change 5Jj in the 
momentum Jj due to the presence of the mode is written as 5Jj/V ^ hq{\uj\'^ — 
|vjp), where q is the wave number of the unstable mode and the particle number 
Nj = JdV\Wj\^ is conserved. Since the total momentum Ji + J2 must be conserved 
in this system, we require the condition 5 Ji + 5 J2 = and 5Jj j^O to exchange 
momentum between the two condensates. This condition is reduced to — 
|vi p = — I W2 P + |v2p 7^ 0, which is satisfied for the unstable modes of the counter- 
superflow instability. Thus the momentum exchange is accelerated as the unstable 
modes are amplified by the counter-superflow instability. 



Critical relative velocity- 
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Fig. 2 Unstable region of counter- superflow. The horizontal axis is the excited mode P nor- 
malized by h/ ah and the vertical axis is the relative velocity Vr normalized by sound velocity 
c = y^grf/m. The gray zone shows Eq. (5). 

3 The dynamics of counter-superflow in a uniform system 

As described in Sec. 1, QT appears in counter-superflow. We discuss the mecha- 
nism of transition to QT by investigating the dynamics of ID and 2D systems. 

Firstly, we investigate the dynamics of ID systems (Fig. 3). Figure 3 (a) shows 
the time evolution of densities of two components. The horizontal axis is x coordi- 
nate and the vertical axis is time. The left of Fig. 3 (a) shows component 1 flowing 
from left to right and the right shows component 2 flowing oppositely, and the rel- 
ative velocity is bigger than V^riticai- According to Fig. 2, when Vr = 1.39V^riticah 
the modes of \P\ < 1.36 are unstable. The density is initially uniform, but the 
waves of |P| < 1.36 are excited in both components by the dynamical instability. 
Through the nonlinearity the waves grow into dark solitons, which travel without 
decaying. The larger the magnitude of the imaginary frequency is, the less time it 
takes for the solitons to appear. The solitons enjoy the repulsive interaction with 
other solitons. The interaction randomizes the direction of propagation of solitons. 
Two BECs exchange momentum, thus momentum decaying as shown in Fig. 3 (b). 
Eventually, in each component, the number of solitons becomes almost constant 
and the numbers of solitons propagating to the right and to the left become almost 
equal, so soliton turbulence appears 

Secondly, we investigate the dynamics in 2D systems. As ID systems, unstable 
waves are excited and grow into dark solitons. Dark solitons are generally unsta- 
ble in systems except for ID systems, so they decay into vortex pairs by snake 
instability^^. The details will be reported elsewhere. 

A 3D system shows the rich dynamics including the transition to QT as well 
as the soliton dynamics like what occurs in ID and 2D systems. Waves are excited 
and grow into dark solitons. The solitons decay into vortices or vortex rings by 
snake instability in 3D systems. These vortices repeat reconnections to develop to 
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(a) (b) 



Fig. 3 (a) Time evolution of density of both components in a ID system. The horizontal axis is 
X coordinate normalized by the coherence length <^ and the vertical axis is time t normalized by 
T = h/gn^. The parameters are n = 1.0 and |V7?|/\^riticai = 1-39 ait = 0, and 7= 0.5. (b) Time 
evolution of momentum. The horizontal axis is time t normalized by T and the vertical axis is 
momentum J normalized by ^/h. 




(a) (b) (c) (d) 



Fig. 4 Surface of the equal density n = OA of one component in 3D two-component QT. The 
parameters are n = 1.0, | V/?|/| V^nticail = 1.66 at r = 0, and 7= 0.5. (a) r/r = 0. Initial state, 
(b) t/z = 39. Waves are excited by the dynamical instability and grow into dark solitons. (c) 
t/z = 43. Solitons decay into vortices by snake instability, (d) t/z = 51. QT appears. 

QT. Then both components are quantum turbulent, where two species of QT are 
interactive (Fig. 4). 



4 Dynamics of counter-superflow in a 2D trapped system 

We propose how to create two-component QT. We address two species of BECs 
with different internal degrees of freedom in a harmonic trap. Let two BECs dis- 
placed oppositely from the center of the potential by applying a magnetic field 
gradient, then reversing the direction of the field gradient and have two BECs 
collide. QT should appear like a uniform system, but it is difficult to study nu- 
merically in a 3D trapped system. Hence we address 2D trapped systems as the 
first trial. We suppose that two cigar-shaped ^^Rb BECs collide with each other. 
In the initial state BECs are separated by a magnetic field gradient (Fig. 5 (a)). 
Reversing the direction of the field gradient, BECs move oppositely and collide 
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Fig. 5 (Left) Time evolution of the density in a 2D trapped system. Parameters are m = 1.5 x 
10~^^kg, a = 5.00nm, and two compoents have same particle numbers N = \ x 10^. The trap 
frequencies are (jOj^ = 27r x 141, cOy = 27ix 445 and co^ is much smaller than co^ and cOy. The 
coordinates x and y are normalized by the trap length Uh — ^JhlmcHx- The value of density is 
shown as the color bar. (a) tuix = 0.0. Initial state, (b) tcOx = 3.9. Waves are excited, (c) tcOx = 6.0. 
Solitons appear, (d) tcOx = 6.4. Solitons decay into vortex pairs. (Right) The time evolution of 
momentum of two components. 



each other. Waves are excited and grow into solitons like a uniform system (Fig. 
5 (b), (c))^. Eventually, solitons decay into vortex pairs (Fig. 5 (d)). Then BECs 
exchange momentum. Thus we obtain the same dynamics as in a uniform system. 
We expect that QT should appear in a 3D system. 



5 Summary 

We studied the nonlinear dynamics of the counter- superflow instability in two 
miscible BECs. The instability induces the momentum exchange between the two 
condensates and the soliton formulations. Then the solitons decay into quantized 
vortices and vortices grow into the quantum turbulence of two superfluids. We 
found that the similar dynamics can occur in a 2D trapped system. Therefore, we 
expect that QT should appear in a 3D trapped system too. We believe that the 
counter- superflow QT offers a new avenue for study of turbulence. The detail will 
be reported soon elsewhere. 
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